Non-equilibrium dynamics of the Holstein polaron driven by an external electric field 
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This work represents a fundamental study of a Holstein polaron in one dimension driven away from the 
ground state by a constant electric field. Taking fully into account quantum effects we follow the time-evolution 
of the system from its ground state as the constant electric field is switched on at t = 0, until it reaches a 
steady state. At weak electron phonon coupling (EP) the system experiences damped Bloch oscillations (BO) 
characteristic for noninteracting electron band. An analytic expression of the steady state current is proposed 
in terms of weak EP coupling and large electric field. For moderate values of EP coupling the oscillations are 
almost critically damped and the system reaches the steady state after a short time. In the strong coupling limit 
weakly damped BO, consistent with nearly adiabatic evolution within the polaron band, persist up to extremely 
large electric fields. A traveling polaron under the influence of the electric field leaves behind a trail of phonon 
excitations absorbing the excess energy gained from the electric field. The shape of the traveling polaron is 
investigated in details. 
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PACS numbers: 63.20.kd, 72.10.Di, 72.20.Ht 

I. INTRODUCTION 

Research in the field of non-equilibrium dynamics of 
complex quantum systems constitutes a formidable theoreti- 
cal challenge. Many advanced numerical techniques, rang- 
ing from exact-diagonalization', expansion using Chebyshev 
polynomials^, time-dependent density matrix renormalization 
group"*, to non-equilibrium dynamical mean field techniques^ 
have been developed to tackle this complex problem. 

More than forty years ago using a path-integral approach 
Thornber and Feynman'' discovered that an electron in a 
parabolic band, driven by the electric field, acquires a constant 
velocity due to emission of phonons. Later approaches to po- 
laron motion in high electric field used Boltzmann equations^, 
the high-field drift velocity was estimated via phonon-assisted 
hopping between different rungs of Wannier-Stark states us- 
ing rate equations^-^. In ref.'' the one-dimensional Holstein 
polaron problem in strong electric field has been mapped 
on a nonstandard Bethe lattice. It has been realized that 
keeping full quantum coherence between many-body states 
is crucial to obtain finite drift velocity for dispersionless op- 
tical phonons. Extensive research of polaron dynamics has 
been conducted within the semiclassical Su-Shrieffer-Heeger 
model to describe properties of conjugated polymers that may 
be used in a variety of applications like molecular electron- 
ics or light-emitting diodes Polaron formation and its 
influence on transport properties has also been investigated 
in the context of DNA molecules within the semiclassical 
Peyrard-Bishop-Holstein modeU^— and other polaron-like 
models22rJ2,. 

Bloch oscillations (BO) represent a fundamental phe- 
nomenon in quantum mechanics where a charged particle in 
a periodic potential exhibits a periodic motion when exposed 
to an uniform external electric field. Since the electrons in 
solids can dissipate energy due to scattering from inelastic 
degrees of freedom on a time scale usually shorter than a 



typical Bloch time it took a long time until the first ex- 
perimental observation of BO was carried out on semicon- 
ducting superlattices2^— and later in optical potentials^"—. 
Nowadays the concept of BO is frequently present in a vari- 
ety of different fields, for instance atomic Bose-Einstein con- 
densates in optical lattices^^, interacting quantum few-body 
systems-''^^" or organic molecules However, the de- 
scription of damping of BO in dissipative medium remains a 
challenging task. 

By choosing the Holstein Hamiltonian as one of the 
simplest model systems describing the interaction between 
a fermion and phonons, we are able to investigate the 
field-induced acceleration of the polaron, which simultane- 
ously dissipates the energy by inelastic scattering on optical 
phonons while maintaining the full quantum nature of the 
problem. Following the time evolution of the ground state 
when the electric field is switched on at time t = 0, we show 
how the polaron reaches the steady state and consequently de- 
velops a constant non-zero velocity. In particular, we calculate 
the steady-state current vs. voltage characteristics of the Hol- 
stein polaron for different regimes of electron-phonon cou- 
plings. 

We discuss the Holstein model in one spacial dimension 
and give a brief overview of the numerical method in the 
second paragraph. In the third we discuss numerical results. 
Here we give special emphasis on the time evolution from the 
ground state towards the steady state by presenting various 
correlation functions in different EP coupling regimes. We 
compare our results with a simple Landau Zener model and 
follow the time evolution of the polaron as it starts propagat- 
ing after switching on the electric field. As a focal point of this 
work we discuss the dependence of the steady state cuiTent on 
the EP coupling and electric field. In the last paragraph we 
give conclusions. 
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II. MODEL AND NUMERICAL METHOD 

We analyze the one-dimensional Holstein model with a sin- 
gle electron, threaded by a time-dependent flux: 

H = -to^(e'^(*)4c,+i+H.c.) 

i 

+ .g^ni(aj + fli) + wo^al'a^, (1) 

i i 

where c] and al are electron and phonon creation operators 
at site i, respectively, and = cjci is electron density. 
ujQ denotes a dispersionless optical phonon frequency and to 
nearest-neighbor hopping amplitude. The dimensionless EP 
coupling strength is A = g'^ /2touJo. The constant electric 
field F that is switched on at time < = enters the Hamilto- 
nian in Eq. [Uthrough the time-dependent phase 6{t) ~ —Ft 
for t > 0. We measure electric field F in units of [io/eoa] 
where cq is the unit charge and a is the lattice distance. We 
furthermore measure time in units of [/i/io]- Unless otherwise 
specified, from here on we set a = eq = /i = io = 1- 

To solve the time-dependent Hamiltonian for a single elec- 
tron coupled to phonon degrees of freedom we use an im- 
proved numerical method, originally introduced in Ref."*^, that 
led to numerically exact solutions of the polaron ground and 
low-lying excited state properties. The method constructs 
the variational Hilbert space (VHS) starting from the single- 
electron Bloch state cj^|0) with no phonons on an infinite lat- 
tice. The VHS is then generated by applying the off-diagonal 
terms of Hamiltonian ([T|| 

{lC''")} = (^'<i„ + i?f)"'^ct|0), ^2) 

where iJkin and Hg correspond to the first and the second 
term of the Hamiltonian in Eq. [T] respectively. Parameters 
Nh and M determine the size of the VHS. In addition, Nh — 1 
represents the maximum distance between the electron and 
the phonon quanta and Nh * M is the maximum number of 
phonon quanta contained in the Hilbert space. The param- 
eter M > 1 (Ref."*^) ensures good convergence in the strong 
EP coupling regime that contains multiple phonon excitations. 
To reach weak coupling regime, A << 1, we introduce an ad- 
ditional parameter iVphmax limiting the maximum number of 
phonon quanta, which enables construction of VHS as large 
as Nh = 40. 

We first solve the Hamiltonian in Eq.[T]for F — 0, i.e., we 
calculate the (zero temperature) polaron ground state^^*^^^. 
Then we switch on the uniform electric field and start the time 
propagation from the initial state using the time-dependent 
Lanczos technique^**. We manage to find numerically accurate 
results of the model away from equilibrium while maintain- 
ing full quantum description of phonons. Since we are deal- 
ing with a single particle in an infinite system, we compute 
time-dependent average of the current operator j(<) = {I{t)), 
where 

Ht) = i (j2 - H-c. j . (3) 



In the case of a time-independent field F, the time-integral of 
the current is directly related to a change of the total energy 

f j{t')dt' = ^h{t)/F^x{t), (4) 
Jo 

where A/i(t) = {H{t)) - {H{t ^ 0)) and x{t) represents the 
travelled distance^. 



III. NUMERICAL RESULTS 

A. Time evolution from tlie ground state towards the steady 

state 

We first present results obtained near the noninteracting 
limit, i.e. at A = 0.01, where j{t) displays damped BO around 
j{t — )■ oo) > 0, see Fig.[Tla). The period and the initial am- 
plitude of BO at small t are consistent with BO of a free elec- 
tron, denoted with thin dashed line in Fig. \Ila). Damping is 
due to inelastic scattering on phonons that is in turn reflected 
in a monotonic increase of the average phonon number (uph) 
with time, as depicted in Fig. [He). Damping is, however, not 
the most important consequence of inelastic scattering. No- 
tably, j{t) approaches a positive steady state current j for 
t > ts ^ Ub and Gis at = 1/5 and F = 1/2, respectively, 
where ts denotes the Bloch oscillation period ts ~ 2tt/F. 
Note that j « jmax = 2. The dependence of the steady 
state current j on F will be discussed further in the text. The 
steady state current as well emerges as a linear dependence 
of the total energy on time: Ah{t) = Fjt + A/iq, see re- 
sults in Fig.llJe), where with increasing t, Ah{t) approaches 
a straight line. In the steady state we as well observe a linear 
increase of (uph) vs. t. When comparing Ah{t) and (uph) in 
the linear regime we find that Ah{t) = ujod{nph)/dt. This 
equality confirms an intuitive expectation that in the steady 
state the total energy gain is entirely absorbed by the lattice. 

On a more technical side we note that to reach a steady 
state, the Hilbert space used in our calculation must contain 
large enough set of excited states that in turn represent the 
reservoir for the absorption of energy. For this reason, differ- 
ent Hilbert spaces were used, depending on the strength of EP 
coupling and the size of F, see as well Caption of Fig.[T] 

At a larger value of EP coupling, A = 0.2, a somewhat 
different physical picture emerges, as shown in Figs. [nt>,d) 
and (f). The main differences can be summarized as: (i) BO 
become overdamped, (ii) j{t) remains positive at all t, and 
(iii) j{t) reaches a steady state after a short time tg < ts. 
Characteristic for a steady state are linear t dependencies of 
(uph) and Ah{t) in Figs.lTJd) and (f), respectively. Common 
to all cases presented in Fig.[Tl is the emergence of a constant 
steady-state current for t > tg. 

In Fig. [2a) and (b) we present current vs. time in the 
strong coupling regime, i.e. at A = 2.0. At F = 1/10 
(Fig. 12 a)) we observe nearly undamped BO as the polaron 
adiabatically follows the polaron band. Regular oscillations in 
(uph) and Ah{t) in Figs.|2c) and (e) portray polaron averages 
nearly identical to their ground state values at corresponding 
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Figure 1: (Color online) j{t) vs. t/tB for two values of F = 1/5 
and 1/2 for ujq = 1, and two distinct values of A: (a) A — 0.01 and 
(b) A = 0.2. Thin dashed line in (a) represents j{t) for A = 0, thin 
horizontal lines in (a) and (b) indicate steady-state values j; {uph) is 
shown in (c) and (d) for the same set of parameters as in (a) and (b), 
respectively; corresponding averages Ah{t) are displayed in (e) and 
(f). Thin dashed line in (e) represents Ah{t) for A = 0. The accuracy 
of time propagation was checked by comparison of the energy-gain 
sum rule, Eq. |4] Parameters, defining functional generator (Eq. |2j 
were iV^ = 40, M = 1, and iVphmax = 6 for _F = 1/5 and 
Nh = 28, M = 1 and Afphmax = 8 for F = 1/2. In this and 
all subsequent figures we used up to A'^st ~ 15 x lO'' states in the 
Hilbert space and A'^stcp = 2000 time steps within each t b ■ Different 
sizes of VHS were used to check the convergence in the thermody- 
namic limit. Thin straight lines represent t oo extrapolations. 




Figure 2: (Color online) j{t) in (a) and (b), {uph) in (c) and (d), and 
Ah{t) in (e) and (f) vs. t/tB for ljq — 1 and A = 2 and three differ- 
ent values of F as indicated in figures. Inset in (a) shows the polaron 
spectrum (ground state, first excited state energy and the continuum 
vs. the wavevector fc) for luq = 1 and A = 2. Note that there is a dif- 
ferent vertical scale used in (e) and (f). We used Nh ~ A'^phmax = 20 
and M = 1 with iV^t = 3 x 10^ 



F ~ 1, Ah{t) approaches a straight line signaling the onset 
of a steady state. 



Determination of the threshold electric field using the 
Landau-Zener formalism 



wavevectors k = Ft = 2nt /ts- The response of the system 
to external field is nearly elastic, since Ah{t = I * ts) ^ 
for any integer value of I. The average current remains indis- 
tinguishable from zero in the largest time interval tested with 
our calculation, i.e. t < 20tB- 

In order to illuminate this behavior we note that in the 
strong coupling limit a large gap A exists in the polaron exci- 
tation spectrum being of the order of ujq. The low-energy po- 
laron excitation spectrum is presented for = 1 and A = 2 
in the inset of Fig.|2ja) where a gap A ^ 0.64 separates the 
polaron band from the excited polaron ban d"^^'^^ — , located 
just below the continuum denoted by the grey area. At small 
F << A there exist exponentially small probability for a 
nonadiabatic transition from the polaron band to the excited 
polaron band or/and into the continuum. 

For sufficiently large F (see Figs.|2tb,d,f) for F = 1/3 and 
1), BO in j{t) lose periodicity even though remnants of BO 
remain clearly visible, and the time averaged current becomes 
finite (nonzero). Additional frequencies appear in (uph) that 
indicate multiple phonon excitations due to polaron transi- 
tions to excited polaron bands. Moreover, the average value 
of (uph) between successive ts intervals increases. The to- 
tal energy Ah{t) as well increases in time. At large field. 



The observed behavior in the strong coupling regime due 
to a large gap in the spectrum resembles Landau-Zener (LZ) 
transition^^'*'", where the probability for tunneling between 
bands in a two level system. 



_ fvFt A/2 
" l^A/2 -vFt 



is given by 



P 



exp 



. (A/2)^ 
vF 



(5) 



(6) 



where A is the energy gap between the two levels and v is the 
velocity. Using Eq.|6]we estimate the threshold electric field 
Fth using 



"th 



(A/2)^ 



(7) 



Such an estimate has been used to determine the dielectric 
breakdown of the insulating half-filled Hubbard model'-^k 
Applying Eq. [T] to the specific case of A = 2 presented in 
Fig. 12] using A ^ 0.64 and v — jmax ~ 0.1, we obtain 
Fth ^ 1-0. LZ formalism gives roughly the correct order of 
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magnitude of Fth, since a noticeable current appears around 
F = 1/3, as seen in Fig.|2jb) and (f) as well as from steady 
state current, presented in Fig. EJc). One should however be 
mindful when considering the LZ formalism. In the polaron 
case the band structure deviates significantly from the ideal 
LZ model with two hyperbolic bands. In the realistic case, 
multiple transitions occur from the polaron band to a contin- 
uum of excited states composed from a polaron and additional 
phonon degrees of freedom. 
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extending behind the moving polaron. Note that the polaron 
is moving from left to right. In the long time limit, 7 is ex- 
pected to be approximately constant, independent of r and t, 
for r sufficiently negative. The average height of the polaron 
tail 7 is due to energy conservation requirement independent 
of A: 



x{t)F ^ {nph)uja x{t)juJo, 



(9) 



therefore 7 ^ F/luq, compare Figs.[3ld) and (h). Note that 
this relation holds only when the system has reached a steady 
state. The length of the polaron tail is given by the expression 
for the travelled distance x{t) = Ah{t)/F. Alt = 8tB and 
A = 0.01 we obtain from Fig.ITJe) Xt = x{8tB) ^ 14.8, that 
fits well with the length of the phonon tail in Fig.[3jd). 

Rather unexpected is the pronounced increase of phonon 
excitations in the forward direction where 7(r) > up to r < 

~ 5 — 7 for alH > presented in Fig. [3] Since time evo- 
lution starts from the ground state at zero temperature, there 
are no phonon excitations present far ahead from the moving 
electron. A substantial forward tail of phonon excitations is a 
consequence of damped BO. Indeed, rf compares well with 
the Stark localization length, i.e. rf ^ Ls = 4/F = 8. Yet 
another intriguing feature in 7(r) emerges as regular oscilla- 
tions in the polaron tail with a period K = uiq/F = 2, clearly 
seen in the small A = 0.01 limit, see Figs.|3jb) through (d). 
At larger A = 0.2 these oscillations become overdamped. 



Figure 3: (Color online) 7(r) for ij;o = 1, A = 0.01 and F — 1/2 in 
(a) through (d) and A — 0.2 in (e) through (h) computed at different 
times. The electric field is switched on at t = 0. Note a vertical scale 
change between (a) and (b) as well as between (e) and (f). Xt in (d) 
represents the travelled distance, see discussion in the text. 



C. Time evolution of polaron 

In Fig. [3] we follow the time evolution of the polaron to- 
wards the steady state at = 1/2. We compute the average 
number of phonon quanta located at a given distance r from 
the electron 

lir) = (^^'"-t+rO-i+r), (8) 
i 

fulfilling the following sum-rule (upk) = J2r 7('')- At i = 0, 
7(r) displays a pronounced peak at the position of the elec- 
tron, i.e. at r = 0, consistent with the shape of the polaron in 
its fc = ground state. After the electric field is switched on, 
j{r) experiences a compelling time evolution with three out- 
standing characteristics: (i) the overall increase of 7(r) with 
time, (ii) development of pronounced asymmetry of 7(r) with 
respect to the electron position at r = 0, and (iii) increased 
amount of polaron excitations in the forward direction. The 
overall increase of 7(r) is consistent with the absorption of 
energy that is deposited in increasing number of phonon exci- 
tations. The asymmetry is a result of a growing phonon tail. 



D. Steady state current 

The focal point of this work is the calculation of the steady 
state current j and analysis of its dependence on F and A. 
In Fig.|4]we present the current-voltage characteristics, i.e. j 
vs. F for different values of A. Note that the upper limit of 
j is given by the current amplitude jmax = 2 in the nonin- 
teracting system. We have limited our calculations to com- 
mensurate values F = ujq/K with K being integer, with 
two exceptions: (i) large F > ujq where we have chosen 
F = 2a;o, Swo, . . . and (ii) results presented with discon- 
nected triangles in Fig. SJa), with details given in the figure 
caption. In the regime A < 0.1, presented in Fig.|4ja), j de- 
creases with increasing F for F > 0.1. Our method does 
not yield steady state results in the regime F < 0.1 due to 
large Stark localization length Ls = Ato/F. Since j = 
for F = as well as in the opposite limit, when F ^ 00, 
there must exist a global maximum value jmax that depends 
on A. For A = 0.1, jmax ^ 0.82, while for A < 0.1, jmax 
is reached somewhere in the interval < F < 0.1, not ac- 
cessible by the present numerical method. Choosing rational 
or even irrational values of wq/F leads to a decrease of j that 
nevertheless remains non-zero even in the latter case. A sweep 
over continuous values of F would lead to spikes in j located 
at integer values of ujq / F, as consistent with observations in 
previous works^"^'^-. 

To gain further insight into the decrease of j with F, we 
plot in Fig. Htb) j/ a/A vs. 1 / \/F and realize that curves ap- 
proximately collapse onto a straight line. The revealed scaling 
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with 1 / VF is a clear signature that we are deahng with a co- 
herent propagation between Stark states with identical total 
energy that are spaced hy K = cjq/F. This is in contrast 
with the assumption of an incoherent hopping between local- 
ized statesi^ , which would predict a dependence j (x 1 /F. 
In turn our derivation, as presented in the Appendix, leads for 
integer if > 1 as well as for luq < W = 4io to a scaling of 
the maximum steady current 



Jo 



F 



(10) 



While Jq cannot be directly compared to the average current 
j, the functional dependence on A and F is in good agreement 
with scaling in Fig.|4|b) that leads to a ^ 0.89 (fit is repre- 
sented by a dashed line). The expression in Eq.[TO]is valid in 
the small A and large F regime (however, F < loq) where j 
decreases due to decreasing overlap between Stark states. 

The scaling breaks down when with decreasing F, j ap- 
proaches the maximum jmax- In Fig. Hfc) we present results 
for larger A G [0.2, 2.0] which enables us to observe the evo- 
lution of j vs. F as the system evolves from the weak EP cou- 
pling (A < 1) towards the strong EP coupling (A > 1) regime. 
With increasing A, the position of jmax shifts towards larger 
values of F while it decreases in its magnitude. The main 
difference between the weak and strong EP coupling regime 
emerges due to increasing energy gap A in the polaron excita- 
tion spectrum^^^^** that for A >> 1 approaches A ^ ojq- Due 
to large A at large X > 1, j remains zero until F ^ Fth- 

We summarize the overview of numerical results with a 
diagram describing different regimes characterized by dis- 
tinct short-time behaviors (after switching on F), presented 
in Fig. lUd). We distinguish four different regimes: (i) The 
regime of damped free particle BO (DBO) for small values 
of g. (ii) Almost critically damped (CD) regime where steady 
state current is reached in a time shorter than or of the order of 
ts, and the oscillations in the current are still visible, however 
j{t) > for any t > 0. (iii) Polaron BO regime (PBO) where 
system evolves nearly adiabatically. Polaron Bloch oscillates 
within the polaron band and damping is exponentially small 
(numerically undetectable). In PBO average current remains 
zero and total energy remains periodic within numerical ac- 
curacy and up to the largest measured time t < 20<b. (iv) 
Damped polaron BO regime (DPBO) where remnants of PBO 
are seen in j {t) while there exists a measurable average cur- 
rent J > within t < 20i_B. 



IV. CONCLUSIONS 

In summary we list our main results. Using a time depen- 
dent Lanczos method we have followed the time evolution of 
the polaron from its ground state towards the steady state af- 
ter the electric field has been switched on. Different sizes 
of VHS have been used to ensure that presented results are 
valid in the thermodynamic limit. Steady state conditions have 
been reached at intermediate to high electric fields and the cur- 
rent vs. voltage characteristics has been plotted for different 
regimes of EP couplings. By calculating the electron-phonon 
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Figure 4: (Color online) Steady state current j vs. F in the weak 
coupling limit in (a) and in the weak to intermediate coupling regime 
in (c), scaling j/\/A vs. 1 / \/F in the weak coupling limit in (b), and 
diagram, presenting different regimes, as described in the text (d). 
Also in (d) circles with error bars indicate positions in the diagram 
where for a fixed g, a maximum value jmax was reached; isolated 
circles, squares and diamonds indicate values used for Figs[T]and|2] 
Disconnected triangles (seven down and one up) in (a) represent j 
using non-integer values (seven rational and one irrational) of 1 jF, 
i.e. F = 5/9, 2/(1 + VS), 5/8, 6/9, 6/8, 7/9, 7/8 and 8/9 at A = 
0.1. Different values of Nh, M and A^phmax were used to ensure that 
error bars, where not specified, are smaller than sizes of the symbols. 



coiTelation function representing the time evolution of the po- 
laron, we show that the absorbed energy in the steady state 
is deposited as an increasing number of phonon excitations 
arranged as a growing tail behind the moving polaron. 



The damped BO can be observed in the extremely weak EP 
coupling limit. In the former case, period of BO Ib — in/F 
should be less than the relaxation time to / 9^ related with the 
emission of phonons. A large gap in the spectrum in the strong 
coupling regime is responsible for observation of nearly per- 
fect BO arising from the polaron motion within the polaron 
band. The breakdown of this quasiadiabatic regime quali- 
tatively resembles the Landau-Zener transition from the po- 
laron band to higher excited states. Analytical estimate for 
the steady state current on the electric field and EP coupling 
constant at large fields is proposed and numerically tested. 
The unusual steady state current vs. electric field dependence, 
j oc y^X/F, valid at large F and small A, reflects the sig- 
nificance of coherent processes for proper description of po- 
laron motion. In contrast, approaches calculating the steady 
state current relying on probabilities for transitions between 
neighboring Wannier-Stark states, mediated by the EP cou- 
pling, yield j cx X/F2^. 
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Appendix: Coherent propagation between Starlt states 

To analyze the propagation of the polaron an alternative ap- 
proach to driven Hamiltonian, Eq.([T]i, is to study eigenstates 
in a constant external electric field, i.e., the electron Hamilto- 
nian is written as 



(A.l) 



Since eigenstates of are localized Stark states we perform 
the transformation to new orthogonal basis. 



ai = 



E 



Wi-lCi, 



(A.2) 



where wavefunctions Wj (being real) are locaUzed in the inter- 
val— L5/2 < j < Ls /2with Ls ~ 4to/-F, and eigenenergies 
ei = Fl + eg. In order to keep constant energy, the particle 
can propagate along the chain only by emitting (absorbing) 
phonons. The novel unperturbed term 



(A.3) 



connects the average displacement in the Stark basis Ae; = 
FAl to phonon generation ujoANph- 

In the following we consider only the simple commen- 
surate case where uo/F — K is integer (A' > 1, i.e. 
F < uq), where electrons can perform coherent hopping be- 
tween Stark states with Al = K keeping Eq constant by 
emitting (or absorbing) single phonon via the coupling term 
H' = g T^iiO'i + ai)- Also restricting phonon frequencies 
to oja < W — 4io we remain in the regime K < Ls- We 
now construct the basis of possible coherent states having the 
same Eq ^ starting with a bare electron state at / = and 
generating novel states by application of H' 



4lo>, 



10), 



(A.4) 



whereby jm denote location of phonons. The matrix elements 
between subsequent states can be evaluated explicitly by ne- 
glecting multiple occupations of sites (being rare for K 1), 
i.e. ji j2 ■ ■ ■ or equivalently simplifying boson factor for 
multiply occupied sites, i.e.. 



^ 5W'j„x-mKWj„-(m-l)K = Tj^^, (A.5) 

which depends within such an approximation only on j,„ = 
jm — raK. We search now for the eigenstates in such a re- 



stricted space in the form 

31 



(A.6) 



where the energies E are obtained solving the system 

Eh^ = Y.T,,b\\ 

31 

= T,,b,+Y,T.j,hi''\ (A.7) 

32 

pL3l---3m _ rj. ,Jl...i,r.-l , ST^ rp ,3l---3m + l 



With inserting the solutions of the Stark problem into 
Eq. dA.SI l. the branching system, Eq. (IA.7b . can be solved quite 
generally. Here, we are interested only in a qualitative behav- 
ior, hence we use the simplification 



-Ls/2 + K<j< Ls/2, (A.8) 



and Tj = elsewhere, where the phase (—1)'^^ emerges from 
fast varying Stark functions Wj in Eq. dA.Sb . 

Solutions of Eqs.( IA.7l ).( IA.8b can be found by an Ansatz 



3l---3m + l 



'ipK 



1 



^{-iy^..+ibi\-^-, (A.9) 



U31---3 
"m+1 



where Ls — Ls — K. The corresponding eigenenergies are 



E = E, 



2.9 



p ~ -^y Lscos{pK), 



(A. 10) 



which leads to group velocities in the tight-binding form Vp cx 



Wo s\tl{pK) with the maximum 



(A.ll) 



The derivation can be made more rigorous taking into ac- 
count the actual Stark wavefunctions Wj and matrix elements 
Eq. dA.Sb . Still it is not expected to change qualitatively the 
scaling of coherent group velocities Vp with the maximum 
tiQ, Eq. dA.l lb . It should be, however, reminded that we did 
not yet match the actual solution Eq. dA.lOl l with the bound- 
ary condition, as determined, e.g., with the first equation in 
the system Eq. dA.7b . Anyhow, it is expected that an eigen- 
state of stationary Hamiltonian, as in Eq. dA.lb . cannot posess 
a finite steady current (solution being a superposition of ±p 
eigenstates). On the other hand, the driven system and the 
time-dependent model, Eq.Q, clearly can generate the cur- 
rent i{t) and in this sense induce solutions with the steady 
current j cx jo = ^^o following Eq. dA.llb . Evidently, more 
rigorous relation between the eigenstates of the stationary case 
and the driven problem is still desired. 
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